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Abstract 

The small size vortex-antivortex pairs proliferation in type II superconducting film is considered 
for the wide interval of temperatures below T c . The corresponding contribution to free energy is 
calculated. It is shown that these fluctuations give the main contribution to the heat capacity of 
the film both at low temperatures and in the vicinity of transition. 
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The important role of the vortex-antivortex pairs in the 2D phase transitions is well 
known. For instance their presence restores the phenomenon of superconductivity in spite 
of the destruction of the long range order by phase fluctuations if. Let us consider the case 
of dirty superconducting film with the thickness d <C £bcs (0) ~ D/T c , where T> = vplt r /3 
is the diffusion coefficient and r tr is the transport collision time. 

Slightly below the BCS critical temperature vortex-antivortex pairing becomes energeti- 
cally profitable and the BKT transition to the state with the finite stiffness takes place. Let 
us recall that this transition is determined by the large scale vortex-antivortex pairs, which 
binding energy has the logarithmic form. The standard description of the BKT transition 
is based on the corresponding logarithmic expression for free energy 
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where n s2 (T) is the superfluid density and £ (T) is the coherence length, which in the vicinity 
of transition has the form: 

e(T) = ^- 

16T c r 

with t = 1 — T/T c as the reduced temperature. Let us underline that this definition differs 
twice from the standard GL expression. The Eq.(JJ) is valid for large size vortex-antivortex 
pairs, when R 3> £. Namely these pairs determine the thermodynamical and transport 
properties of the 2D superconductor in the vicinity of the BKT transition. 

In this letter we will demonstrate that below the immediate vicinity of the BKT transi- 
tion, in contrast to the current opinion , the essential role in fluctuation properties of 2D 
superconductor play the small size vortex-antivortex pairs. The cornerstone of presented the- 
ory consists in the fact that the energy of the vortex-antivortex pairs turns zero when the 
distance between their centers becomes less than £ . As the consequence such "cheap" 
pairs become "affordable" for thermal fluctuations and their proliferation takes place even 
at so low temperatures when other quasiparticles are frozen. In result namely this small size 
vortex-antivortex pairs give the dominant contribution with respect to other fluctuations in 
the Ginzburg-Landau region below Tbkt and even prevail over the BCS exponential tail in 
the heat capacity of the 2D superconducting system at low temperatures. 

In order to take into account the specifics of the mentioned fluctuation processes let us 
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start from the general expression for the partition function 



Z = J®A(r) J £A(r)*exp 



F (A (r) , A (r) 
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with F (A (r) , A (r)*) as the Ginzburg-Landau functional, which is valid in the vicinity of 
T c , and then we will generalize the result obtained for the case of low temperatures. In con- 
trast to the usual Ginzburg-Landau type long wave-length approximation the calculation 
of the functional integral here has to take into account also the vast variety of the order 
parameter functions A p (r) corresponding to the specific realization of the vortex-antivortex 
pairs in the film. In other words side by side with the usual Ginzburg-Landau order pa- 
rameter fluctuations we will take into account additionally some specific short wave-length 
fluctuations. 

Let us first separate the partition function Z of the superconducting film without fluc- 
tuations: 

Z = Z ■ Z(fiy 

The partition function Z^ we will calculate in the gas approximation. Namely we will 
assume that the main contribution to it appears from the small pairs neglecting their overlap. 
Hence 

z (/0 - z p 



Z, = / dA s (r) / dA 5 (r)- exp I , (2) 



where 

T 

is the contribution to the partition function of the isolated single pairs of all possible sizes 
5 ■ £, with < (J < 1. The power S/ (vr£ 2 ) in Eq.(j2J) takes into account the combinatorial 
factor corresponding to the independent formation of such pairs. The order parameter 
As (r) must have at the distance 2<5£ two zeros of the opposite vorticity (i.e. the total 
vorticity calculated at large distances is equal to zero). 

As the next step we disregard the axial asymmetry of the vortex-anivortex pair. We will 
see below that the characteristic size of the pairs which mainly contribute to the partition 
function r e ff £ what will justify our gas approximation (5 < 1. The corresponding to 
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Eq.(J2J) free energy functional is 



F P = »dJ dh [-r| A 5 (r) | 2 + g|ft_A, (r) | 2 + A s (r) 



(3) 



Here u is the density of states and <9_ = d/dr—2ieA, ( (x) is the Riemann zeta-function. 
One can demonstrate [3] that the energy of As (r) profile with two zeros at the fixed distance 
2<5£ is majorized by that one of the disc of radius 5C, and below we will estimate the latter. 

The GL functional minimization procedure is equivalent to the solution of the 2D 
Ginzburg-Landau equation with the boundary condition As (r = 8£) = 0. Let us assume 
the order parameter in the form 

As (p) = A (T) / (p) 

and write down the dimensionless GL equation: 

1 d ( d \ „ 1 1 



Here p — r/£ is the dimensionless radius, while A (T) is the BCS gap. 

Let us investigate the solutions of Eq. (j3J) in two limiting cases: small p ~ 5 1, and 
large p — > oo. In the region of small p ~ 8 1 the function / (p) — > and it can be found 
as the solution of the linearized Eq. (J3J): 

/ (p) = d Jo (p/W) + C 2 iV (p/v^) , (5) 

where Jo and A^o are the Bessel and Neumann functions correspondingly. 

When p > 1, the solution of Eq.(jH) can be looked for in the form / (p) = 1 — fx (p) . 
Substituting this expression to Eq.(jU) one finds 

f(p) = l-Mp) = l-CsK (p), (6) 

where K is the Bessel function. The constant term corresponds to the usual long wave- 
length order parameter contribution which will appear in the result of our consideration side 
by side with the contributions of the small vortex-antivortex pairs. 

One can see that the asymptotic of Eq.fjHJ), being calculated for p > 1 nevertheless extends 
to the region of small values of p. It is why we can match both solutions Eq. (|SJ) and Eq. (jHJ) 
at p ~ 8 <C 1. In order to do this let us use the expressions for all Bessel functions J , N and 
K at small arguments j4|. The comparison of the coefficients at logarithmic term and of 



constants gives two equations for coefficients Ci,2,3 the third equation is determined by the 
boundary condition / (5) = which corresponds to our disc model for the vortex-antivortex 
pair. Finally 



Ci = Sit + 
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Here 7 = exp (C) , where C — 0, 577215... is the Euler constant. 

Now we can find the energy of disc (vortex-antivortex pair) as the function of its dimen- 
sionless radius S. In order to do this one can use the general expression (JHJ) and the fact that 
A = A (T) / satisfies the nonlinear Eq. P|). Using this and performing the integration of 
the derivative by parts one can find the energy of vortex-antivortex pair: 



F P (6,T) 
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Let us stress that namelly presence of the large logarithm in denominator of Eq. ((7|) does 
such localized vortex-antivortex type fluctuations of the order parameter more energetically 
favorable than the almost homogeneous (long wave-length) Ginzburg-Landau ones. 

Now let us perform the functional integration in Eq.(j2J). Vortex-antivortex pairs prolif- 
erate below T c in fluctuation way and in our model this is equivalent to the account of the 
discs of the different radii: 



Z„ = 2 Sd8 exp 



F P (S,T) 
T 
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The integral in Eq. (JHJ) may be carried out by the method of the steepest descent what gives 

4-2 1 /4 vr V 2 (B(T)^ 1/4 
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and finally, the corresponding contribution to the free energy is 

TS 



F P {T) = -TlnZ^t 



ne (T) 



2^T) 1 m 
T 4 T 



(10) 



5 



This expression is valid in the assumption B (T) ^> T which is necessary for the applicability 
of the steepest descent method. Looking on Eq. (jHJ) and recalling the definition of the 
characteristic for superconducting film 2D Ginzburg-Levanyuk number 

Gi « = 2 ^W!Z « 1 <»> 

one can see that this requirement is equivalent to r > Gi^d), i-e. our consideration is valid 

for temperatures below the critical region of transition. 

At temperatures close to transition, in the GL region 1 ^> r ^> Gi^d) the main temper- 
ature dependence of Eq. (jTUJl originates from functions £ (T) and A (T) . Corresponding 
vortex-antivortex pairs contribution to heat capacity is 



C p (T -> T c ) = 



t>ir/ tr . Or 2 




t , 2n 4 udVr 
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One can see that the differentiation of the second term in the square brackets of the Eq. 
(II 2j) gives nothing else as the well known contribution to heat capacity occurs due to the 
GL long wave length order parameter fluctuations |5j]. 

<> (1 » r » = ^g„T c f . (13) 



Nevertheless the main fluctuation contribution to the heat capacity of a 2D superconducting 
film (d <C £ (t)) results from the first term of the Eq. (fT2"j) . corresponding to the protagonists 
of this letter, i.e. short wave length vortex-antivortex type order parameter fluctuations: 



C\ fl) >(l^T^Gi {2d) ) = -^-yvT c ^— — j . (14) 



One can see that the contribution (fPfj) dominates over ()13j) in the entire Ginzburg-Landau 
region r > Gi^d)- 

In the critical region r < Gi^d) the main contribution to heat capacity belongs to the 
classical BKT large vortex-antivortex pairs jsj]: 

(r « Gi m ) ~ -vTJSd (^) 3/2 

which matches with Eq. (fPfj) at r ~ Gi^d) ■ Let us note that the negative sign of the 
vortex-antivortex fluctuation corrections means that the mean field heat capacity jump 
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overestimates its true value and the vortex-antivortex pairs fluctuation contribution smears 
it out. 

Now let us turn to discussion of the region of low temperatures. Unfortunately we cannot 
write down the exact functional for the superconductor's free energy here, but it is almost 
evident that the main changes in it with respect to the Ginzburg-Landau one at T ~ T c 
occur through the different temperature dependence of £ (T) and A (T) at low temperatures. 
Hopefully the region of validity of the Eq. (fTUj) in this way can be considerably extended 
and, at least as the estimate, we will still use it just substituting the GL parameters £ (T) 
and Ao (T) by their exact BCS values. 

The microscopic BCS equation for A (T) is well known 



ln^ = 27rTV (— 

Wn>0 \ 
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with uo n = 2irT(n + 1/2). In the limiting cases 
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In the region of low temperatures T<T C the temperature corrections to the order parameter 
are exponentially small by the parameter exp (— A (0) /T) . 

The coherence length £ (T) in the most general case can be defined as the pole of the 
linear response operator kernel for the corresponding correction to modulus of A |7( . In the 
simplest case of dirty superconductor the corresponding equation takes form 



E 
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. K + CO] (v^ + AfUT) - X?/ [2t? (T)] 
In the limiting cases of low and high temperatures 
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where 6 = 0.76595 is the solution of the equation 



T^T C 
T < T c 
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Now we are ready to calculate the contribution of the vortex-antivortex fluctuations to 
the heat capacity C^ v ~ a ^ (T <C T c )of the2D superconducting film (d <C £ (0)) at low temper- 
atures. The values £ (T) and A (T) here can be assumed to be temperature independent: 
£ (T) « £ (0) = 2fl ^ (0) and A (T) « A (0) = vrT c /7 and differentiation of the Eq. (JTUJ) 
results in 

(-(«-«) ( T « T c ) 7ttC (3)H ffjl^ Q g7 / Gi(2^\ 1/2 ^ (i5) 



AC (T c ) 8(Z tr p F ) 7 2 V T d ' V * 

where £ = T/T c <C 1 and 

7<(3) 

is the BCS jump of heat capacity at T c . 

Contribution (fTH|) requires a special discussion. First of all let us stress its power depen- 
dence on temperature and unusual growth with the decrease of temperature. Comparison of 
Eq. (j!5|) with the BCS expression for heat capacity of superconductor shows that already 
at T < T c j In GiT^s the former exceeds the latter, i.e. nonmonotoneous temperature de- 
pendence of a heat capacity can be expected. Naturally the formal divergence of Eq. (JTSJl 
at very low temperatures has not be taken seriously: at low temperatures the quasiclassical 
approximation used for evaluation of the vortex-antivortex energy breaks down and 

these excitations probably acquire some small gap (~ A 2 /Ep) in the spectrum (in complete 
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analogy with the excitations in the vortex core in 2D superconductors 

In conclusion let us summarize the results obtained. We have demonstrated that the 
possibility of small vortex-antivotrex pairs proliferation in 2D superconducting film results 
in appearance of the specific contribution to its heat capacity. Being gapless such excitations 
result in the appearance of the non-exponential heat capacity of superconducting film at 
low temperatures which turns to dominate the BCS term. At high temperatures the small 
vortex-antivortex pairs contribution to heat capacity dominates over all other fluctuation 
corrections in the entire GL region below T c . Finally, only within the critical region the 
large BKT type vortex-antivortex pairs become of the first importance and determine the 
thermodynamics of the 2D superconducting film. 
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